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I.  Introduction 

The  absorption  and  backscatter  of  light  in  the  underdense  blowoff 
of  a  laser  produced  plasma,  as  well  as  a  limitation  of  electron  thermal 
energy  flux  all  play  important  roles  in  laser  fusion.  This  paper 
constitutes  a  one  dimensional  fluid  simulation  of  these  and  other 

relevant  processes.  It  is  a  continuation  of  our  earlier  work  in 
1  2 

this  area  *  but  it  models  the  dynamics  of  the  interaction  rather 
than  the  steady  state.  Also  included  in  this  work  is  a  study  of  the 
effects  of  Brillouin  backscatter.  Although  our  work  applies  in 
principle  to  any  laser  produced  plasma,  it  is  particularly  motivated 

3_6 

by  experiments  at  the  Naval  Research  Laboratory.  Experiments 

there  have  demonstrated  that  backscatter  through  the  lens  (apparently 

Brillouin  backscatter)  can  be  a  very  important  process  for  structured 

3-4 

laser  pulses  at  sufficiently  high  intensity.  Also,  more  recent 
experiments  with  longer  pulses^  (pulse  duration  of  about  3  nsec)  have 
shown  very  high  absorption  even  at  intensities  where  inverse  brems- 
strahlung  is  not  effective.  For  instance  at  an  irradiance  of 
7  x  10^  W/cm2,  the  fractional  absorption  is  about  60%.  Also  experi¬ 
ments  at  NRL^  and  elsewhere'7  have  shown  that  absorption  is  enhanced 
when  the  target  is  in  the  focal  plane  of  the  lens.  It  has  always 

been  our  opinion  that  resonant  absorption  alone  cannot  account  for 

2 

all  of  the  absorption.  Reference  2  reviews  absorption  measurements 
done  at  many  laboratories  and  comes  to  this  conclusion.  The  addi¬ 
tional  mechanism  which  we  propose  is  ion  acoustic  turbulence  driven 

by  the  return  current.  That  mechanism  takes  place  in  either  an 
Note:  Manuscript  submitted  July  20,  1979. 
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unmagnetized  or  magnetized  plasma.  The  fluid  simulations  reported 
here  attempt  to  model  this  process.  A  very  much  oversimplified 
picture  is  that  ion  acoustic  turbulence  in  the  underdense  plasma 
begins  to  take  over  where  classical  inverse  breamsstrahlung  becomes 
ineffective.  For  instance,  the  solid  line  in  Fig.  1  is  a  schematic 
of  the  total  electron-ion  momentum  exchange  collision  frequency  as  a 
function  of  laser  irradiance  for  a  Nd  laser  produced  plasma.  The 
dotted  line  is  a  plot  of  classical  inverse  bremsstrahlung,  while  the 
dashed  line  is  a  plot  of  the  effect  of  ion  acoustic  instability. 
Roughly  speaking,  ion  acoustic  turbulence  at  high  temperature  gives 
perhaps  half  as  much  absorption  as  inverse  bremsstahlung  at  low 
temperature. 

One  of  the  principal  conclusions  of  this  work  is  that  the  absorp¬ 
tion  physics  falls  into  two  different  regimes.  In  short  single 
pulse  experiments,  resonant  absorption  is  the  dominant  process,  but 
it  gets  a  strong  boost  from  the  ion  acoustic  turbulence.  Typically 
the  total  absorption  is  about  40%  with  resonant  absorption  accounting 
for  25%  and  ion  acoustic  turbulence  for  15%.  Most  of  the  light  is 
specularly  reflected,  although  stimulated  Brillouin  backscatter  plays 
a  role  also.  On  the  other  hand,  for  long  pulses  or  structured  pulses, 
the  physics  is  usually  dominated  by  a  relatively  complicated  inter¬ 
play  between  backscatter  and  ion  acoustic  turbulence.  Both  of  these 
processes  depend  upon  the  quasi-linear  and  nonlinear  behavior  of  a 
plasma  instability,  so  the  physics  is  necessarily  more  speculative 
than  for  instance  absorption  by  inverse  bremsstrahlung.  This  modeling 
is  more  complicated  in  a  fluid  simulation  but  can  nevertheless  bring 
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a  more  complete  picture  of  absorption  in  laser-plasma  interactions. 

Our  treatment  of  these  processes  is  internally  self  consistent  and 

is  based  upon  well  documented  theories,  simulations  and  experiments. 

We  will  now  briefly  contrast  our  approach  to  the  physics  of  the 

underdense  plasma  with  particle  simulations  and  other  fluid  simulations. 

Particle  simulations  have  been  particularly  useful  in  studying 
8—10 

resonant  absorption  , strongly  driven  backscatter  instabilities  in 
a  laser  produced  plasma^  and  current  driven  ion  acoustic  instabil¬ 
ities  in  infinite  homogeneous  plasmas.  The  problem  with  the  latter 
two  processes  is  that  they  are  slow  compared  to  the  laser  frequency 
so  that  particle  simulations  get  very  expensive.  For  instance  the 
simulation  of  Brillouin  backscatter  in  Ref. 11  ran  to  t  =  7  x  103/fi 
(35,000  time  steps)  with  V  /V  ^7.  However,  this  is  only  3.5  psec 
for  a  Nd  laser  produced  plasma,  and  the  total  length  corresponds  to 
only  10  ym,  meaning  a  density  gradient  scale  length  of  1-2  pm. 

Clearly,  a  particle  simulation  for  more  realistic  times,  lengths  and 
smaller  VQg/Ve  fields  would  be  astronomically  expensive.  Our  approach 
sacrifices  the  detailed  description  of  the  instabilities.  Its  effect 
is  now  modeled  by  anomalous  transport  in  a  fluid  code.  However  we 
gain  tremendously  in  the  range  of  parameters  that  can  be  studied. 

For  instance,  for  a  70  psec  Nd  laser  pulse,  we  model  the  entire 
underdense  plasma  for  the  entire  pulse  duration  with  .  i  than  two 
minutes  computing  time  with  a  TI-ASC. 

The  other  approach,  used  to  study  laser-plasma  interactions,  is 
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fluid  simulations  (characterized  by  LASNEX  ).  Our  simulations 
have  more  in  common  with  them  than  they  do  with  particle  simulations 
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but  there  are  still  Important  differences.  Specifically  our  treatment 
of  flux  limitation  is  quite  different.  In  our  model,  a  large  thermal 
flux  excites  an  instability.  This  instability  tends  to  limit  the  flux, 
and  gives  rise  to  several  other  anomalous  transport  effects.  All  of 
these  transport  effects  can  be  self  consistently  related  to  each  other 
through  the  turbulence  level.  To  estimate  this  level  we  appeal  to 
previous  experiment,  simulation  and  theory.  Thus  flux  limits  are 
derived  in  this  work,  not  specified.  Also  this  work  incorporates  a 
model  for  Brillouin  backscatter.  Also,  as  discussed  in  the  next 
section,  our  approach  to  electron  thermal  energy  transport  is  quite 
different  from  that  used  in  LASNEX.'*'^ 

The  outline  for  the  remainder  of  this  paper  is  as  follows.  In 
the  next  section  we  discuss  the  anomalous  physical  processes  in  our 
simulation  models.  Section  III  summarizes  the  fluid  equations  which  we 
solve  numerically.  Section  IV  gives  results  of  test  runs  on  the  code. 
Section  V  presents  a  number  of  different  simulations  and  compares  with 
experiment.  Specifically  we  look  into  absorption  for  short  pulses, 
long  pulses,  double  structured  pulses  and  absorption  as  a  function  of 
Z  and  laser  wavelength.  Comparisons  with  recent  NRL  experiments  are 
shown.  Section  VI  suggests  future  experiments  and  Section  VII 
summarizes  our  conclusions. 

II.  Anomalous  Physical  Processes 

A.  Ion  Acoustic  Instability  Generated  by  Return  Current 


The  problem  of  ion  acoustic  instability  in  a  plasma  which 


conducts  energy  has  been  discussed  in  a  series  of  papers  which  deal 
with  both  unmagnetized^’^  and  magnetized^ ’ ^  plasmas.  In  this  work 
we  confine  our  study  to  unmagnetized  plasmas.  The  basic  idea  is  that 
if  the  electrons  conduct  energy  but  carry  no  net  current,  the  flux  of 
energetic  particles  in  the  direction  of  Q  (the  heat  flux)  must  be 
balanced  by  a  return  current  of  low  velocity  particles  in  the  opposite 
direction.  This  return  current  can  excite  ion  acoustic  waves  which 
also  propagate  opposite  to  Q.  There  are  four  principal  physical 
effects  of  this  instability.  First  of  all  it  reduces  the  energy  flux, 
or  in  other  words,  reduces  the  electron  thermal  conduction.  Second, 
it  gives  rise  to  an  anomalous  energy  exchange  between  electrons  and 
ions,  with  the  ions  gaining  energy  and  the  electrons  losing.  Third, 
it  gives  rises  to  an  electric  field  which  generates  the  return  current 
(this  latter  process  is  not  directly  calculated  in  this  work  because 
we  deal  only  with  the  total  momentum,  not  with  the  momentum  of  each 
species  separately) .  These  three  processes  are  discussed  in  detail  in 
Ref.  16.  From  the  quasi-linear  theory  of  a  one  dimensional  spectrum  of 
ion  acoustic  turbulence  which  is  antiparallel  to  Q,  the  calculation  of 
thermal  conduction  and  energy  exchange  is  straightforward  and  the 
results  are  given  in  the  next  section. 

Fourthly,  the  ion  density  fluctuations  give  rise  to  enhanced 

1  2  16  18 

absorption  of  the  laser  light.  *  ’  ’  Making  a  Lorentzian  approxi¬ 

mation  to  the  electron  distribution  function,  the  anomalous  collision 
frequency  for  absorption  of  the  laser  light  is  from  Ref.  16 


where  E  is  the  electric  field  of  the  laser  light,  ft  is  its  frequency, 
and  $(k)  is  the  fluctuating  potential  due  to  the  ion  acoustic  wave  at 
wave  number  k.  If  the  laser  light  incident  normal  to  the  plane  of  the 
target  and  the  ion  acoustic  spectrum  is  one  dimensional,  k*E  =  0.  Of 
course  the  turbulent  spectrum  is  not  one  dimensional  but  will  exist 
within  a  cone  of  angle  0  with  respect  to  the  return  current.  Since 
our  results  do  depend  upon  ecJj/T^  and  8,  it  is  important  to  pick  these 

parameters  as  accurately  as  possible.  Fortunately  we  can  draw  on  a 

19-26  27-33 

great  deal  of  experimental  data  for,  and  particle  simulations 

of  current  driven  ion  acoustic  instabilities.  Also  there  is  a  great 

34**38 

deal  of  recent  theory  which  examines  ion  trapping,  resonance 

broadening  or  anomalous  transport  and  which  basically  supports  these 
data.  Our  basic  assumption  is  that  an  ion  acoustic  instability  driven 
by  a  return  current  will  have  the  same  fluctuation  level  and  angular 
spread  as  one  driven  by  a  real  current. 

19 

We  will  now  review  some  of  these  data.  Mase  and  Tsukishima 

examined  a  current  driven  ion  acoustic  instability  and  diagnosed  it  by 

microwave  scattering.  They  confirmed  that  the  waves  obeyed  the  linear 

dispersion  relation  and  found  5^"  <v  0.1.  Also,  they  measured  the 

e 

angular  spread,  about  the  direction  of  current,  of  the  wave  numbers 
of  the  fully  developed  turbulence.  For  -30°  <  0  <  30°,  (as  far  as  they 
measure  in  angle)  they  find  a  spectrum  which  is  nearly  uniform  as  a 
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function  of  angle.  Finally  they  confirm  that  the  anomalous  electron 

collision  frequency  is  proportional  to  which  supports  the  quasi- 

e 

linear  theory  we  invoke. 

Also,  there  are  measurements  of  ion  acoustic  turbulence  in 

collisionless  shocks.  The  turbulence  at  k  -  kp  itself  was  first 

20 

identified  by  CO^  laser  scattering  in  the  shock  front.  Later  experi- 
21 

ments  measured  the  angular  width  of  this  spectrum  and  found  6  ^  50° 

22 

to  60°.  A  more  recent  experiment  produced  ion  acoustic  turbulence 

by  driving  a  current  through  a  collisionless  plasma  column.  The  linear 

dispersion  relation  was  verified,  and  it  was  found  that  for  a  strongly 

driven  system  (VQ  ^  0.15  Vg)  the  fluctuation  level  was  about  15%.  The 

authors  of  Ref.  22  also  compared  the  measured  collision  frequency  with 

38  23 

that  predicted  by  Choi  and  Horton.  Hamburger  also  showed  that  ion 

acoustic  instability  could  be  driven  in  a  toroidal  system  with  high 

current.  However  no  measurements  of  fluctuation  spectrum  or  angular 

24 

spectrum  are  given.  Even  in  a  weakly  driven  system  ,  ion  acoustic 

24 

turbulence  can  still  be  excited.  In  the  experiment  of  Slusher  et  al  , 


the  current  velocity  in  He  at  Tg  =  4eV  exceeds  the  threshold  velocity 

by  a  factor  of  two  or  less.  Also  impurity  ions  (protons  in  this  case) 

were  found  to  have  a  strong  stabilizing  effect.  Nevertheless  ion 

acoustic  waves  with  k  'v  k^/2  were  observed  in  accord  with  linear 

theory.  By  CO2  laser  scattering,  they  found  'v  0.03  and  an  angular 

e 

width  of  about  30°.  Finally,  experiments  at  University  of  California 
at  Los  Angeles  also  found  well  diagnosed  ion  acoustic  waves  with 
k  'v  kp/2  and  which  grow  to  a  fluctuation  level  of  ^  0.1.  These 
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ion  acoustic  waves  could  be  driven  either  by  an  electron  current  or 

2  6 

by  the  decay  of  an  electron  plasma  wave. 

We  now  turn  to  a  discussion  of  particle  simulations  of  current 

27-29 

driven  ion  acoustic  turbulence.  Lampe  et  al  have  shown  with  a 

one  dimensional  particle  simulation  that  if  a  current  is  perpendicular 

to  an  applied  weak  magnetic  field,  ion  acoustic  fluctuations  grow  as  in 

an  unmagnetized  plasma.  These  simulations  also  confirmed  that 

k  ^  k^/2.  The  instability  was  found  to  saturate  by  ion  trapping  (see 

next  subsection)  with  0.1  <  ^  <  0.2.  Also  these  simulations  confirmed 

e 

that  the  electron  dynamics  were  described  by  quasi-linear  theory. 

There  are  also  several  simulations  of  current  driven  ion 

30-32 

acoustic  turbulence  in  two  dimensions  .  In  Ref.  30,  resonant 

quasi-linear  heating  was  again  confirmed  for  the  electrons.  Once  the 

system  saturated  by  ion  trapping,  with  0.1  <  <  0.2,  the  turbulence 

e 

had  a  wide  angular  spectrum,  8  ^  55°  -  60°.  Reference  31  claimed  ion 

trapping  is  not  responsible  for  saturation.  However  they  see  ion 

acoustic  turbulence  with  ^  0.1,  kAp  *  1  and  a  spectrum  of  angular 

e 

width  0  •v  30°.  Also  they  confirm  quasi-linear  behavior  for  electrons. 
32 

Biskamp  et  al  also  confirmed  that  in  a  two  dimensional  particle-in¬ 
cell  simulation,  a  current  drives  ion  acoustic  waves  unstable.  They 

found  a  broad  angular  spectrum  with  ^  >  0.1  and  kAQ  <  1,  and  the 

e 

electrons  described  by  quasi-linear  theory.  Another  similar  simula- 

33 

tion,  but  in  three  dimensions  confirms  these  results. 

We  have  modeled  the  two  dimensional  nature  of  the  spectrum  by 
assuming  that  (k)  is  constant  within  a  cone  of  angle  0  about  the 
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direction  of  the  return  current,  and  zero  outside.  Also  we  model  the 

spectrum  with  a  single  k.  A  choice  0  =  k  ^  kp/2  and  ^-  =  0.1  seems 

e 

to  be  reasonably  consistent  with  the  data  which  we  have  cited. 

Let  us  finally  note  that  if  the  plasma  is  magnetized,  the  ion 
wave  has  its  wave  vector  perpendicular  to  both  B  and  VT^,  so  k  can  be 
parallel  to  the  electric  field  over  part  of  the  laser  spot.  Thus  not 
only  can  the  absorption  be  enhanced,  but  also  there  is  no  need  to  make 
any  assumption  concerning  angular  width  of  the  spectrum.  Our  study, 
which  neglects  B,  then  gives  rise  to  minimum  anomalous  absorption. 

B.  Brillouin  Backscatter 

We  now  discuss  the  Brillouin  backscatter  instability  in  an 

39  40  41 

inhomogeneous  plasma.  ’  ’  To  start,  we  briefly  review  the  results 

of  the  linear  theory  of  Brillouin  backscatter  in  an  inhomogeneous 

plasma.  Since  the  flow  velocity  in  the  underdense  plasma  is  generally 

supersonic,  we  assume  that  the  group  velocity  of  the  reflected  wave 

and  acoustic  wave  are  in  the  same  (outward)  direction.  However,  even 

if  the  group  velocities  have  opposite  sign,  the  result  of  the  theory  is 

39  42 

not  different.  ’  In  the  steady-state,  the  linear  equations  for  the 
incident  and  reflected  wave  are 


V 
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where  F  denotes  the  square  root  of  the  action  density,  a  subscript  r(s) 
denotes  the  reflected  (sound)  wave  and  V  is  the  appropriate  group 


velocity.  The  quantity  1q  is  the  growth  rate  for  the  Brillouin  back- 

scatter  instability  in  the  absence  of  damping  of  the  reflected  and 

,  40 

sound  wave. 


where 


(3a) 


(3b) 


where  E.^  is  the  amplitude  of  the  electric  field  of  the  incident  laser 
light.  Finally  K  is  the  reciprocal  of  the  square  of  the  phase  mismatch 
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where  k  is  the  wave  vector  of  the  appropriate  wave  (for  Brillouin 

backscatter  k  =  k^  and  k  ~  2kj),  and  v  is  the  Landau  damping  rate 
r  i  s  i  s 

of  the  sound  wave  which  will  be  specified  in  the  next  section. 

For  the  case  v  =  0,  the  reflected  wave  travels  back  towards 
s 
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the  laser,  and  is  spatially  amplified  for 

W"  Lc  5  Vk(V  V  W2  (4) 

gs  gi 

and  the  total  energy  amplification  is 

exp  t2iTY02/K(VgsVgr)  !/2] .  (5) 
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On  the  other  hand,  if  the  sound  wave  is  strongly  damped,  it  is 
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amplified  in  a  region  of  approximate  length 


1*1  <  Lc  =  kT  (6) 

s 

and  the  total  energy  amplification  is  still  given  by  Eq.  (5). 

Let  us  now  discuss  how  to  model  this  process  in  our  fluid 
code.  First  imagine  that  the  grid  size  Ax  is  larger  than  2  Lc«  Then 
a  reflected  wave,  which  is  in  phase  here,  is  amplified  only  within 
this  particular  cell  and  not  in  neighboring  cells.  Thus  the  cells  act 
independently  of  each  other.  If  we  know  the  thermal  level  for  the 
reflected  wave,  we  can  figure  out  how  much  it  is  amplified  in  each 
cell.  This  amount  is  then  subtracted  from  the  incident  laser  light 
and  2Vg/c  times  this  power  goes  into  ion  heating  (the  actual 
algorithm  will  appear  in  the  next  section).  Clearly,  if  Ax  >  2Lc 
the  total  backscattered  light  is  independent  of  grid  size,  since  if 
the  grid  is  twice  as  large,  the  thermal  level  of  reflected  light  in 
the  grid  cell  is  twice  as  great  and  so  is  the  total  reflection  in  that 
grid  cell.  As  an  initial  reflected  noise  level,  we  usually  assume 
that  for  the  grid  size  chosen,  ten  energy  e-folds  in  a  cell  gives 
total  reflection  there. 

Now  we  turn  to  the  physics  of  the  more  complicated  case  of 
Ax  <  2L^.  In  this  case,  the  reflected  wave  amplifies  coherently  from 
one  cell  to  the  next.  To  treat  this  would  require  keeping  track,  in 
each  cell,  of  many  different  reflected  waves  which  are  resonant  in 
many  different  cells.  Because  this  was  not  the  usual  case,  we  treated 
the  cells  as  independent  of  each  other.  Our  approach  is  to  first 
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assume  an  energy  amplification  in  the  cell  given  by  Eq.  5  and  then 


reduce  the  total  amplification  by  Ax/2Lc.  This  gives  the  right  answer 
if  the  amplification  is  coherent  for  2Lc/Ax  cells,  but  of  course  this 
algorithm  is  certainly  not  as  valid  as  that  for  the  case  of  Ax>2Lc» 
Actually,  our  calculations  usually  show  that  Ax  >21^. 

As  we  will  see,  if  this  linear  theory  is  used,  the  result  is 
that  Brillouin  backscatter  is  an  extraordinarily  strong  process  and  it 
dominates  the  dynamics  (that  is,  the  laser  light  is  almost  totally 
reflected)  for  nearly  all  gradient  scale  lengths  and  incident  laser 
light  intensity.  Particle  simulations  have  also  confirmed  that  while 
in  the  linear  regime,  Brillouin  backscatter  can  play  a  dominant  role 
even  if  the  gradient  scale  length  is  of  order  of  a  free  space  wave- 
length^Thus  a  crucial  problem  is  to  see  whether  any  process  reduces 
the  strength  of  Brillouin  backscatter  instability. 

We  consider  a  nonlinear  reduction  in  the  growth  rate  due  to 
ion  trapping  by  the  ion  acoustic  wave.  To  start,  we  will  show  that 
as  the  ion  wave  grows,  it  is  forced  into  a  nonlinear  regime.  For 
temporal  growth,  where  vg  =  0  each  photon  of  incident  wave  lost 
produces  one  photon  of  reflected  light  and  one  phonon  of  ion  acoustic 
wave.  For  spatial  growth,  the  situation  is  different  because  the 


group  velocity  of  the  photon  is  much  greater  than  the  group  velocity 

42  44 

of  ion  acoustic  wave.  ’  Here  an  incident  photon  decays  into  one 


reflected  photon  and  V  /V  reflected  phonons.  From  this,  it  is  a 

gr  gs 

straightforward  matter  to  calculate  the  saturated  value  (that  is 
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assuming  total  reflection) 


13 


a 


given 


where  Z  is  the  ion  charge.  Clearly  the  saturated  value 


(?) 

^  ej sat 


in  Eq.  (7a)  is  much  larger  than  that  required  for  ion  trapping  for  any 
reasonable  value  of 


V  /V  .  Also,  (lr^\  given  by  Eq.  (7b)  is  compared 

08  /e<t>\  V  e/Sat 

to  the  maximum  value  of  (Tp-]  given  by  Eq.  (8a)  and  (8b).  It  is  always 

\  e/tr 

larger  by  at  least  a  factor  of  5  or  10,  and  often  much  more.  There¬ 
fore,  before  the  photon  can  be  strongly  backscattered,  the  ion  acoustic 
wave  is  driven  into  a  strongly  nonlinear  regime. 

To  describe  the  backscatter  in  this  regime,  we  invoke  recent 
37  45 

theoretical  work.  ’  This  work  does  not  utilize  an  enhanced  damping, 
but  rather  a  modification  to  the  growth  rate  in  the  absence  of  damping. 
The  idea  is  that  there  is  some  power  input  to  the  ion  acoustic  wave. 

In  the  linear  regime,  this  power  input  causes  the  ion  acoustic  wave  to 
grow  at  the  proper  linear  growth  rate.  However  once  the  ions  begin  to 
trap,  a  small  increase  in  amplitude  of  the  ion  acoustic  wave  traps  a 
few  more  ions.  Since  ions  abruptly  gain  a  tremendous  amount  of  energy 
when  they  are  trapped,  the  power  input  goes  principally  into  acceler¬ 
ating  ions  and  only  slightly  into  causing  wave  growth.  However  the 
wave  does  not  stop  growing  completely  because  the  power  input  is  still 
there.  Nevertheless,  the  growth  rate  does  abruptly  drop  once  ions 
begin  to  trap.  The  calculation  of  the  growth  rate  is  given  in  Ref.  37 
for  a  wave  which  grows  in  time  in  an  homogeneous  plasma  and  in  Ref.  45 
for  a  wave  which  grows  in  space  in  an  homogeneous  but  bounded  plasma. 

In  either  case,  it  was  found  that  this  process  could  be  described  in 
terms  of  a  large  reduction,  by  a  factor  of  between  five  and  ten,  of  the 

growth  rate  y (that  is  the  growth  rate  if  vg  =  0).  We  invoke  these 
results.  If  a  comparison  of  Eqs.  (7)  and  (8)  shows  that  trapping  is 
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important,  we  reduce  yq  and  also  Lc>  according  to  Eq.  (3c)  by  this 

large  factor.  Depending  on  whether  this  new  reduced  is  greater  than 

or  less  than  half  a  grid  cell,  we  treat  this  backscatter  as  incoherent 

or  coherent  from  cell  to  cell  as  discussed  previously. 

To  summarize,  we  first  check  to  see  whether  the  Brillouin 

backscatter  instability  is  dominated  by  ion  trapping.  Where  it  is  (and 

it  always  is  in  the  simulations  we  have  done)  we  adopt  previous  theo- 
37  45 

ries  ’  which  show  that  the  growth  rate  yq  is  substantially  reduced. 
Let  us  note  that  the  approach  used  here  is  very  different  from  that 
used  in  Ref.  43.  Here  we  rely  on  reducing  yq»  and  also  Lc  according  to 
Eq.  (4).  The  latter  invokes  an  increase  of  (according  to  Eq.  (6)) 
until  it  is  longer  than  the  plasma  size.  In  a  future  work,  we  will 
give  a  more  detailed  discussion  of  these  two  approaches. 

C.  The  Problem  of  Energetic  Electrons 

Other  numerical  studies  of  this  problem,  principally  LASNEX 
have  assumed  that  the  electron  energy  is  transported  mostly  by  ener¬ 
getic  electrons  which  are  treated  separately.  The  basic  idea  is  that 

the  classical  mean  free  path  is  much  longer  than  any  characteristic 

21  —3 

scale  length.  For  instance  if  Tg  =  lOkeV  and  n  =  10  cm  ,  the 
electron  mean  free  path  is  roughly  one  millimeter.  Thus  a  classical 
fluid  model  is  not  valid  for  these  electrons.  However  separate  trans¬ 
port  of  these  electrons  is  a  very  complicated  procedure  and  to  our 
knowledge,  there  is  no  generally  agreed  upon  procedure  for  it. 

The  approach  that  we  adopt  is  very  different.  We  use  a  single 
fluid  model  and  justify  it  by  a  large  anomalous  collision  frequency. 

For  instance,  if  the  return  current  excites  an  ion  acoustic  instability 
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6(f) 

with  fluctuation  amplitude  ,  the  anomalous  collision  frequency  is 

e  6  ed>  1 

given  roughly  by  w  (e<*r  )2.  Therefore,  if  e/£  ^  10  the  electron 

pe  e  ie 

21  -3 

mean  free  path  for  a  10  keV  plasma  at  n  =  10  cm  is  roughly  2 
microns.  Thus  a  single  fluid  approximation  is  now  much  more  reason¬ 
able.  Even  if  the  absorption  process  gives  rise  to  nonthermal  electron 

distributions  (for  instance  resonant  absorption  almost  certainly 
8—11 

does  ),  if  the  transport  of  these  energetic  electrons  is  inhibited, 
a  single  fluid  theory  is  probably  still  valid  everywhere  except  right 
at  the  position  where  the  nonthermal  electrons  are  accelerated.  At 
these  points  the  thermal  temperature  we  calculate  will  be  some  average 
of  the  super thermal  and  thermal  temperature  actually  produced. 

Thus  our  philosophy  is  to  explore  the  consequences  of  a  single 
fluid  description  of  the  electrons  and  see  where  it  leads.  For  one 
thing  we  find  that  the  temperature  in  the  underdense  plasma  is  gener¬ 
ally  much  greater  than  the  thermal  electron  temperature  predicted  by 
LASNEX.  Our  thermal  temperature  in  the  underdense  region,  in  fact 
agrees  much  better  with  both  the  superthermal  electron  energy  of  the 
LASNEX  simulations  and  shown  by  experiment.  Of  course  there  will 
still  be  runaway  tails  on  the  electron  distribution  function,  but  now 
they  perturb  rather  than  dominate  the  electron  transport.  We  find 
that  this  model  does  give  reasonable  results.  Another  great  advantage 
is  that  it  is  a  very  straightforward  theory  and  is  very  economical 
to  study. 
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III.  Fluid  Equations 


As  mentioned  in  our  earlier  work  ,  the  advantages  of  a  time- 
46 

dependent  approach  lie  not  only  in  the  information  it  yields  on  the 
dynamics  of  the  absorption  processes  but  also  in  the  choice  of  boundary 
conditions  away  from  the  interaction  region  (which  is  itself  part  of 
the  solution).  In  the  present  work,  boundary  conditions  were  chosen 
both  in  or  near  the  solid  and  in  the  underdense  region,  at  densities 
several  orders  of  magnitude  above  or  below  the  critical  density.  The 
fluid  equations  given  now  describe  the  evolution  of  mass  density, 
total  momentum,  electron  and  ion  temperature 
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where  n  is  the  electron  density,  m  the  electron  mass,  M  the  ion  mass, 
v  the  fluid  velocity,  Tg  and  T^  the  electron  and  ion  temperature,  Z 
the  average  ion  charge,  Qg  and  the  electron  and  ion  thermal  flux, 
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—  L  4  .  —  •  .m, 


v>c  and  v  the  classical  and  total  electron-ion  collision  frequency 

(v  =  v  +  v  where  the  second  subscript  stands  for  anomalous),  cm 
c  an  Te 

the  anomalous  electron-ion  energy  exchange  term,  E^  and  the 

incident  and  reflected  laser  field  intensities,  P  the  resonant 

res 

absorption  heating  power,  g  the  backscatter  heating  term.  P^rem 
is  the  bremsstrahlung  radiation  loss  term  (at  the  temperature 
achieved  in  the  runs,  the  elements  under  consideration  are  completely 

ionized).  The  expression  for  these  latter  terms  are: 

1/2 
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v  n  Z  In  A 

c  = 


3.5  105  T  3/2 
e 


(16) 


where  vg  is  the  ion  sound  speed  and  A  will  be  defined  later.  The 
parameter  a  will  be  discussed  also  in  the  next  section;  x£  is  the 
location  of  the  critical  surface.  It  is  important  to  note  that  in 
this  model  all  the  absorption  processes  heat  up  the  electron  thermal 
distribution  (because  of  the  fluid  modeling)  .  There  are  no  non- 
thermal  tail,  nor  a  double-temperature  for  the  electrons  in  this 
work.  It  has  been  shown  that  ion-acoustic  turbulence,  which  is  an 
important  part  of  the  absorption  physics,  does  not  lead  to  very  flat 
tail  but  seems  to  deposit  its  energy  for  electrons  with 
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Of  course  this  is  also  true  of  classical  inverse 


v  =  2-3V  . 

e 

bremsstrahlung. 

The  electron  thermal  flux  term  is  defined  as: 

9T 


Qe  =  -  n 


dx 


with 


K^-1  =  K  -1  +  K, 
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(17) 
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an  ce 

whereas  the  ion  thermal  flux  term  remains  classical 

3T . 

Qi  ■  *Kci »  ir  <19) 

The  equations  for  the  incident  and  reflected  electric  fields  read: 

1/2 
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where  A  is  the  Brillouin  backscatter  term.  It  is  defined  as 
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(21) 


where  Ax  is  the  local  grid  size,  YQ/6  the  nonlinear  spatial  growth 

rate.  The  factor  of  6  models  the  effect  of  ion  trapping  as  discussed 

in  the  previous  section  and  the  factor  Max  (Ax,  L£)  accounts  for  the 

fact  that  L  (the  scattering  length)  may  be  longer  than  a  grid  cell, 
c 

In  addition  to  solving  Eq.  (20),  we  also  solve  the  equation  for  the 
reflected  wave  with  A  set  equal  to  zero.  This  gives  the  specularly 
reflected  light  intensity.  Substracting  the  specularly  reflected 
light  from  the  total  reflection  gives  the  Brillouin  backscatter. 
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An  important  step  in  the  process  of  solving  these  equations  is 
the  determination  of  the  unstable  regions  both  for  the  ion-acoustic 
and  Brillouin  backscatter  instabilities.  Unlike  the  classical  terms, 
the  anomalous  and  backscatter  terms  are  non-zero  only  when  some 
criteria  for  these  instabilities  have  been  met.  We  turn  now  to  these 
criteria  and  discuss  also  how  they  are  implemented  in  the  numerical 
solution. 

A.  Ion-acoustic  Instability 

1  2 

In  our  steady-state  work,  ’  we  solved  for  the  amplitude  of  the 
unstable  waves  and  we  got  the  fluctuation  density  level  as  part  of 
the  solution.  Because  of  the  extra  physics  included  in  the  present 
work  and  because  of  the  added  constraints  due  to  the  time-dependent 
approach,  a  simplified  treatment  for  the  determination  of  the  effects 
of  the  ion-acoustic  turbulence  was  adopted.  The  level  of  fluctuation 
was  chosen  to  be  a  certain  prescribed  value  (0.1  in  all  the  runs 
except  in  the  test  runs  reported  in  the  next  section) .  This  level 


~  ’v  0.1  is  supported  by  numerous  experiments  and  particle  simu- 
e 

lations  of  current  driven  acoustic  turbulence  as  discussed  in  the 


(23) 


(24) 


where  0  is  the  angle  between  the  unstable  ion  wave  cone  and  the  laser 
light  direction.  It  was  discussed  in  the  previous  section  and  set  to 
tt/3  in  our  calculations.  Also  k  is  the  average  unstable  ion  wave 
vector  and  taken  to  be  k^  /2. 

The  growth  rate  of  the  ion  acoustic  wave  is  given  by 
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The  procedure  for  including  the  anomalous  transport  coefficients  due 
to  ion-acoustic  fluctuations  is  the  following:  for  every  point  in 
space  x,  starting  with  the  old  value  for  Qg  (defined  at  the  previous 
time  step,  whether  or  not  it  corresponded  to  an  unstable  situation), 
we  calculate  y.  We  then  calculate  a  new  value  for  Qg,  based  on  the 
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value  obtained  for  y.  We  then  recalculate  y  to  check  that  its  value 


is  consistent  with  the  calculation  of  Q^.  Iterations  are  performed 
until  consistency  is  achieved.  We  define  also  a  flux  limitation  factor 

Qe 

f  =  — rr-r  which  will  be  shown  in  the  results  section.  This  flux 
nmv  i 
e 

inhibition  factor  is  obtained  as  part  of  the  solution  and  as  such,  is 
space  and  time  dependent.  At  low  density,  f  is  generally  limited  by 
the  instability,  while  at  high  density  the  plasma  is  usually  stable 
with  the  classical  value  of  f  being  sufficiently  small.  Thus  the  flux 
limit  f  is  not  arbitrarily  chosen,  but  is  calculated  self  consistently. 
If  the  f  calculated  this  way  is  too  large,  an  over-riding  flux  limit  of 
f  =  1.5  is  invoked.  However,  this  upper  flux  limit  is  almost  never 
used;  at  most  for  a  few  points  in  space  for  only  the  first  few  time 
steps  of  the  calculation. 

B.  Stimulated  Brillouin  Backscatter  Instability 

The  same  type  of  approach  has  been  adopted  for  the  stimulated 
Brillouin  backscatter  instability.  The  Brillouin  backscatter  term  A 
is  set  to  zero  everywhere  at  the  beginning  of  each  time  step  and 
differs  from  zero  only  when  the  following  double  criterion  is  met 
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where  vr^sj»  ^gr(s)  are  t^ie  damping  rate  and  group  velocity  of  the 
reflected  (sound)  wave 


We  have  consistently  checked  that  the  ion  trapping  conditions  (8a)  and 
(8b)are  satisfied.  We  have  also  checked  that  vs<wL  so  that  we  expect 
Eq. (2b) to  be  valid. 

The  factor  of  10  appearing  in  Eq.  (21)  for  A  was  introduced 

to  provide  physically  a  smooth  transition  between  the  cases  of 

complete  reflection  and  no  reflection  at  all.  It  simply  means  that 

10  growth  lengths  are  needed  for  complete  reflection  in  a  cell  if 

the  cell  size  Ax  is  larger  than  2L  .  Unlike  the  previous  instability, 

c 

no  iteration  on  time  was  necessary  but  we  note  that  y0  and  A  depends 

on  the  value  of  the  incident  field  (through  V  ).  The  incident  field 

os 

is  affected  in  turn  by  the  Brillouin  backscatter.  The  strong  spatial 
coupling  between  Brillouin  backscatter  and  the  laser  light  must  be 
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reflected  in  the  algorithm  and  as  shown  in  Appendix  A,  it  was  found  to 
be  crucial  for  an  acceptable  solution  to  this  problem. 

IV.  Test  Runs 

In  the  previous  section,  the  equations  were  given  along  with  the 

procedures  to  calculate  various  other  relevant  physical  quantities 

which  complete  the  description  of  our  model.  It  turns  out  that  three 

quantities  needed  to  be  specified  initially;  the  effects  of  these 

quantities  on  the  solution  is  now  investigated.  These  three  quantities 

are  the  density  fluctuation  level  the  fractional  absorption  due 

e 

to  resonant  absorption  for  any  ray  reaching  the  critical  density,  and 

the  non-linear  saturation  of  the  Brillouin  backscatter  light. 

A.  Level  of  Ion  Acoustic  Density  Fluctuations 

1  2 

In  our  previous  steady-state  studies  *  (unmagnetized  and 

magnetized),  we  solved  for  the  unstable  ion  acoustic  node  amplitudes 

and  as  a  result,  we  obtained  calculated  values  for  The  charac- 

^e 

teristics  of  the  solution  generated  were  the  following:  the  maximum 
total  fluctuation  amplitude  reached  a  level  of  about  0.1;  the  fluctua¬ 
tion  level  maximum  spanned  a  region  of  a  few  tens  of  microns  for  the 
unmagnetized  case  (smaller  for  the  magnetized  case  but  we  are  not  con¬ 
cerned  with  this  case  in  our  study);  the  same  modes  made  up  for  most 
of  the  total  fluctuation  level  except  for  some  cases  where 
different  modes  came  into  play  at  different  locations. 

From  these  results,  we  took  0.1  as  our  standard  value  for  the 
density  fluctuation  level.  In  this  section  only,  we  varied  e<k/T 

e 

between  0.05  and  0.2.  Results  have  been  obtained  for  a  short  pulse 
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and  a  double-structured  pulse.  These  types  of  initial  conditions  will 

be  described  in  detail  later.  The  peak  laser  irradiance  in  the  former 

case  is  10^  W/cm2  and  for  the  latter  is  3  x  10^  W/cm2.  We  find  that 

for  both  cases,  the  total  absorption  increases  with  ^  as  well  as 

Ae 

the  maximum  electron  temperature  achieved  during  the  course  of  the 

run.  For  the  10^  W/cm2  case,  total  absorption  varies  from  31%  for 

e$/T  =  0.05  to  63%  for  e^/T  =  0.2  (42%  for  e<J>/T  =  0.1).  The 
e  e  e 

maximum  electron  temperature  ranges  from  3.4  keV  to  11  keV  with  5  keV 

for  the  intermediate  case.  A  very  similar  pattern  emerges  for  the 
15  ,  n 

3  x  10  W/cm^  case.  The  differences  are  even  more  accentuated  since 
the  total  absorption  goes  from  19%  to  69%  while  the  maximum  electron 
temperature  varies  from  4.8  keV  to  19  keV.  The  conclusions  for  this 
test  are  thus  the  following: 

-  The  results  do  depend  on  the  value  chosen  for  the  ion 

acoustic  fluctuation  level  ^  . 

Ae 

-  Because  of  its  dependence  on  this  parameter,  this  study 
points  out  the  importance  of  this  absorption  mechanism  and  its 
definite  role  in  laser  light  interaction  with  targets,  not  only  for 
the  enhanced  absorption  but  for  the  anomalous  transport  and  flux 


inhibition. 

-  Because  of  the  dependence  of  results  on  it  is 

e 

important  that  this  parameter  be  chosen  as  accurately  as  possible.  A 
value  'v  0.1  is  in  reasonable  agreement  with  our  earlier  steady 

state  calculations.  Also  it  is  consistent  with  many  simulations. 


experiments  and  theories  of  current  driven  ion  acoustic  instability. 
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B.  Fractional  Resonant  Absorption 

Resonant  absorption  for  optimum  angle  of  incidence  has  been 
shown^’^  to  be  around  50%  for  p-polarization.  Averaging  over  angles 
of  incidence  and  over  polarizations,  a  reasonable  number  for  absorp¬ 
tion  due  to  this  process  would  be  around  20-30%.  In  the  remainder  of 
this  work,  we  have  chosen  30%  (a  =  0.3)  to  be  the  fraction  of  the 
laser  light  which  is  absorbed  when  it  reaches  the  critical  density. 

In  this  section,  we  varied  o  between  0.2  and  0.4.  Results  have  been 
obtained  for  the  same  kind  of  pulses  and  initial  conditions  used  for 
the  previous  tests.  Without  going  into  any  detail,  the  total  absorp¬ 
tion  changes  less  in  these  tests  than  in  the  previous  ones.  This  is 
due  mainly  to  the  smaller  range  of  parameter  variations  in  this  case 
[factor  of  2  variation  vs.  a  factor  of  16  variation  for  (e<|>/Te)2]. 

The  increase  in  total  absorption  follows  closely  the  increase  in 
resonant  absorption  and  in  both  cases,  the  backscattering  remains 
almost  unchanged.  This  indicates  that  for  each  case,  the  same  regime 
applies.  The  increase  in  maximum  electron  temperature  is  very  minimal 
between  the  20%  and  40%  resonant  absorption  case,  of  the  order  of 
1  keV.  In  conclusion,  because  the  range  of  uncertainty  in  the  reso¬ 
nant  absorption  level  is  not  as  large  as  for  the  ion  turbulence,  the 
results  depend  only  weakly  on  the  exact  level  of  resonant  absorption 
used  in  the  calculations.  For  this  reason,  it  is  our  feeling  that  its 
modeling  is  fairly  accurate  and  that  the  results  shown  later  properly 
take  into  account  this  absorption  mechanism  (or  any  other  mechanism 
occurring  near  the  critical  surface). 
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C.  Non-linear  Saturation  of  Brillouin  Backscatter 

In  the  previous  sections,  non-linear  saturation  of  the 

Brillouin  backscatter  was  discussed  at  length  since  it  constitutes  one 

novelty  of  this  work.  In  this  section,  we  report  on  two  different 

kinds  of  test.  First,  a  run  was  made  for  the  case  of  a  hydrogen 

plasma  without  reduction  in  growth  rate  (due  to  ion  trapping)  but  with 

43 

the  ion  heating  due  to  backscatter  instability.  The  initial  density 
gradient  scale  length  was  taken  to  be  lOp  (the  steeper  the  gradient  is, 
the  less  backscatter  is  expected  and  observed  experimentally)  and  the 
laser  irradiance  was  10  W/cm  .  The  backscatter  reaches  quite  high 
values  for  this  case  and  is  equal  to  48%  after  70  psec.  The  total 
absorption  is  18%  which  is  quite  low  compared  to  similar  runs  with 
similar  laser  irradiances  as  will  be  seen  in  the  results  section.  Since 
ion  heating  gives  rise  to  maximum  reduction  in  growth  rate  for  Z  =  1, 
there  will  be  even  more  backscatter  for  higher  Z  (most  runs  have 
Z  =  2.6  corresponding  to  CH?  targets  in  the  result  section).  This  test 
shows  clearly  the  importance  of  introducing  ion  trapping  in  our 
Brillouin  backscatter  model.  Another  series  of  test  involved  changing 
the  reduction  factor  in  the  linear  growth  rate  in  order  to  mock  up  the 
non-linear  one.  A  factor  of  6  is  commonly  used  in  all  production  runs. 
This  factor  was  changed  to  5  and  7  in  the  following  runs.  For  a  laser 
irradiance  of  3  x  10  W/cm  and  a  double-structured  pulse  (the  defini¬ 
tion  of  this  term  will  be  given  later),  the  total  absorption  varied 
from  28%  for  the  reduction  factor  of  5  to  58%  for  the  reduction  factor 
of  7.  The  corresponding  backscatter  was  found  to  be  respectively 
equal  to  55%  and  17%.  The  difference  between  the  5  and  6  reduction 


factor  cases  was  quite  small.  The  conclusion  for  these  tests  is  then 
that  ion  trapping  is  very  important  in  order  to  simulate  results  frot 
actual  experiments  and  that  the  level  of  backscatter  computed  is 
qualitative  only. 

D.  Grid  Size 

A  standard  test  on  grid  size  (which  controls  the  time  step 

through  the  Courant-Friedrichs-Levy  condition  for  an  explicit  scheme) 

was  performed  in  order  to  check  the  validity  of  the  calculations. 

Besides  its  numerical  significance,  this  test  -  where  the  grid  size 

was  halved  -  is  quite  important  because  the  physics  near  the  critical 

density  can  be  quite  sensitive  to  the  grid  size.  For  example,  it  has 

14 

been  pointed  out  by  Rosen  et  al  that  they  could  not  see  some  predict¬ 
able  physical  effects  near  the  critical  surface  because  of  their  too 
low  resolution  in  that  region.  Since  in  most  of  our  calculations,  the 
smallest  grid  size  is  just  below  2y  near  the  critical  surface,  it  was 
quite  interesting  to  see  changes  -  if  any  -  brought  by  a  ly  cell  size 
in  that  region.  Changes  observed  were  very  small,  at  most  of  the  order 
of  a  percent  or  so  on  most  quantities.  This  is  in  fact  what  was 
expected  because  the  Nd  laser  light  wavelength  is  ly  in  vacuum  only, 
and  increases  near  cut-off  due  to  WKB  enhancement.  Features  due  to 
the  ponderomotive  force  like  the  lower  density  shelf  as  reported  in 
other  works^’ was  well  observed  as  well  as  the  upper  density 
shelf  in  some  cases.  Emphasis  will  not  be  placed  on  these  features  in 
the  next  section  since  this  is  not  the  main  subject  of  our  study. 
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V.  Production  Runs 


Runs  were  made  for  various  laser  pulse  intensities,  for  various 
initial  density  gradients,  for  various  elements  and  for  several  laser 
wavelengths.  For  these  runs,  we  report  on  total  absorption,  resonant 
absorption,  backscatter,  specular  reflection,  maximum  electron  temper¬ 
ature,  flux  limitation  f  =  Q/nmVe3  at  the  critical  density  and  maximum 
ion  expansion  velocity  between  0.1  n^  and  nc  during  the  course  of  the 
runs.  The  laser  pulse  intensities  vary  in  peak  value  (maximum  inten¬ 
sity)  and  in  temporal  distribution.  In  the  following,  short  pulses 
will  refer  to  a  Gaussian  pulse  shape  with  a  FWHM  of  75  psec  for  Nd  and 
750  psec  for  light;  a  long  pulse  essentially  looks  the  same  until 

it  reaches  its  peak  (50  psec  for  Nd;  500  psec  for  CC^)  when  the 
intensity  is  kept  a  constant.  The  initial  density  gradients  belong  to 
two  different  categories:  a)  the  constant  density  gradient  below 
critical  density;  b)  the  double  structure  density  gradient.  The  idea 
behind  this  classification  is  the  following.  At  t  =  0,  it  is  assumed 
that  the  surface  material  of  the  target  has  already  been  ablated  by  the 
laser  light  for  a  very  short  time  compared  to  the  pulse  length.  In 
case  of  the  short  pulse,  it  is  assumed  that  a  lOp  density  gradient 
scale  length  for  Nd  laser  light  is  present  when  the  pulse  is  turned  on 
(lOOu  for  CO2  light) ;  for  a  long  pulse  it  is  assumed  that  a  lOOp 
density  gradient  scale  length  had  time  to  form  in  front  of  the  target 
(lOOOp  for  CO2  laser).  These  define  the  constant  density  gradient 
cases.  For  the  double  structure  density  gradient,  it  is  assumed  that 
the  target  has  been  hit  by  a  prepulse  sometime  before  the  arrival  of 


the  pulse  under  consideration  (main  pulse).  In  this  particular 

situation,  we  choose  typically  a  30p  density  gradient  scale  length 

below  n  down  to  0.1  n  ,  and  a  100p  density  gradient  scale  length 
c  c 

below  0.1  n  .  This  density  structure  is  in  reasonable  agreement  with 
c 

that  measured  by  interferometry  for  a  target  irradiated  with  a  pre- 
3 

pulse.  No  such  runs  were  performed  for  CO^  laser  lights. 

This  double-structure  density  gradient  is  an  artifice  for  simu¬ 
lating  a  prepulse  because  it  proved  to  be  inefficient  to  run  a  real 
prepulse  case  with  the  pulse  and  prepulse  separated  by  2  or  3  nsec. 
Also  2D  effects  might  become  important  on  this  time  scale.  For  all 
the  runs  in  this  section,  the  following  parameters  were  set  at  the 
values  of  0.3  for  a  (fractional  resonant  absorption  once  the  light 
reaches  critical  density),  and  0.1  for  e<)>/Te.  Also  all  the  runs 
reported,  except  when  mentioned  otherwise,  have  been  made  for  poly¬ 
ethylene  targets  with  a  maximum  Z  of  2.67. 

A.  Single  Short  Pulses 

In  these  runs,  because  of  the  shorter  density  gradient  scale 

length,  backscatter  is  expected  to  be  unimportant.  This  is  what  the 

experiments  show,  in  any  case.  A  typical  solution  obtained  at  the 

15  2 

peak  of  the  laser  pulse  for  3  x  10  W/cm  is  shown  in  Fig.  2.  On  the 
electron  density  curve,  let  us  note  a  lower  density  shelf  forming 
above  0.5  nc«  The  electron  temperature  profile  shows  a  maximum  in  the 
lower  density  region  around  7  keV.  The  ion  temperature  presents  some 
structure.  The  maximum  near  the  critical  density  is  due  to  anomalous 


energy  transfer  between  electrons  and  ions.  The  bumps  at  lower 
densities  are  usually  a  residual  of  backscatter  in  that  region.  Most 
of  the  plasma  is  unstable  to  ion-acoustic  instability  and  the  flux 
inhibition  factor  shows  a  maximum  in  the  underdense  region  (minimum 

inhibition).  Above  0.1  n  ,  it  is  first  negative  which  corresponds  to 

3Te  -3 

a  negative  value  of  Qg  (-j^—  >  0) .  It  then  goes  as  low  as  10  ,  goes 

up  again  as  the  density  increases  and  reaches  a  maximum  value  above 

critical  density  and  tapers  off  around  0.3.  This  constitutes  a  fairly 

typical  behavior  for  f.  Notice  that  the  ion  acoustic  turbulence  not 

only  limits  the  thermal  flux  in  the  underdense  plasma,  but  also  limits 

it  up  to  about  lOy  into  the  overdense  plasma  also.  At  sufficiently 

high  density,  the  thermal  flux  becomes  classical  again.  The 

number  for  f  which  is  reported  in  all  the  following  figures  is 

taken  at  n  =  n  .  Finally,  the  fluid  velocity  (which  is  not  shown) 
c 

O 

reaches  a  maximum  value  of  1.4  10  cm/sec(corresponding  to  50  keV  ions) 

at  n  =  0. 1  n  .  At  this  point  in  time,  the  absorption  for  this  case 

was  38%,  25%  of  which  was  due  to  resonant  absorption.  A  fraction  of 

9%  had  been  backscattered  whereas  52%  had  been  specularly  reflected. 

The  variation  of  these  quantities  with  laser  intensities  is  shown  in 

Fig.  3.  From  this  figure,  we  see  that  the  total  absorption  remains 

roughly  constant,  around  40%,  over  2  orders  of  magnitude  in  intensity. 

The  resonant  absorption  is  responsible  for  more  than  half  of  the  total 

16  2 

absorption  (around  25%).  Only  above  10  W/cm  ,  does  the  backscatter 
increase,  the  resonant  absorption  drops  leading  to  a  decrease  in  total 
absorption.  This  calculated  absorption  of  40%  is  in  reasonable 
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agreement  with  short  pulse  experiments  on  low  Z  targets  at  both  NRL^9 
12 

and  Livermore.  Figure  4  shows  maximum  electron  temperature,  ion 
blowoff  velocity  and  flux  limit  at  n  =  n^  as  a  function  of  irradiance. 
Also  plotted  as  individual  points  on  this  figure  are  values  of  hot 
electron  temperature  which  were  measured  at  the  Naval  Research 
Laboratory"^  for  75  psec  Gaussian  pulses  on  polyethylene  targets. 

B.  Long  Pulses 

Our  simulations  of  long  pulse  experiments  in  Nd  laser  produced 

plasma  generally  show  that  there  is  an  initial  transient  during  which 

the  plasma  settles  down.  This  treatment  typically  lasts  for  tens  of 

picoseconds,  and  is  usually  characterized  by  strong  Brillouin  back- 

scatter.  The  duration  of  this  transient  period  increases  with  irra- 

1 A  2 

diance.  For  instance  it  is  less  than  50  psec  for  10  W/cm  and  is 

15  2 

about  70  psec  for  10  W/cm  .  After  this,  the  plasma  settles  down 
into  a  much  steadier  state  (for  typically  50-150  psec  until  the  end  of 
the  run)  as  far  as  absorption  is  concerned.  The  ions  however  accel¬ 
erate  throughout  the  run,  so  maximum  ion  velocity  may  actually  be 
larger  than  shown  here.  All  of  the  results  we  quote  here  are  for  the 
latter  part  of  the  pulse. 

The  total  absorption  is  larger  in  these  long  pulse  experiments 
than  in  the  case  of  short  pulse  experiments,  probably  because  the  ion- 
acoustic  turbulence  has  much  more  space  to  grow  on  these  large  density 
scale  lengths.  Also,  because  the  absorption  takes  place  in  this 
manner,  less  light  actually  gets  to  the  critical  region  and  the 
percentage  for  resonant  absorption  is  smaller  than  for  single  short 
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pulse  experiments.  The  absorption  physics  is  then  dominated  by  the 
complicated  interplay,  in  the  underdense  region,  between  classical 
and/or  anomalous  absorption,  and  backscatter. 

In  Fig.  5  are  plotted  absorption,  resonant  absorption,  back¬ 
scatter  and  specular  reflection  as  a  function  of  irradiance  from 
12  15  2 

3  x  10  -  1.5  x  10  W/cm  .  Experimental  results  on  total  absorption 

5 

are  shown  as  circles  for  the  NRL  long  pulse  experiment  ,  and  as  squares 

for  the  Sandia  long  pulse  experiment. At  low  irradiance,  the 

absorption  is  classical  and  no  instabilities  are  excited.  At  higher 

irradiance,  both  ion  acoustic  turbulence  and  backscatter  are  important. 

As  a  rule  of  thumb,  the  transition  from  classical  to  anomalous 

1 A  2 

behavior  occurs  above  10  W/cm  ,  at  which  point  the  total  absorption 

is  80%,  but  absorption  due  to  ion  acoustic  turbulence  is  13%  (which 

was  measured  by  redoing  the  calculation  with  return  current  driven  ion 

acoustic  instability  turned  off),  resonant  absorption  is  16%  and 

14  2 

fraction  backscattered  is  3%.  At  7  x  10  W/cm  ,  the  anomalous 

absorption  (measured  in  the  same  way)  is  54%  out  of  a  total  of  66% 

with  12%  resonant  absorption  and  15%  backscatter. 

In  Fig.  6  are  shown  maximum  electron  temperature,  ion  blow- 

off  velocity  and  flux  limit  as  a  function  of  irradiance.  Note  that  at 
14  2 

7  x  10  W/cm  ,  the  maximum  temperature  7.5  keV  compares  reasonably 
well  with  the  nonthermal  tail  temperature  of  10  keV  quoted  in  Ref.  5. 
The  picture  which  emerges  from  these  calculations  is  that  at  high 
irradiance,  the  nonthermal  tail  on  the  integrated  x-ray  distribution 
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comes  from  a  thermal  plasma  at  low  (roughly  critical)  density,  while 
the  lower  temperature  x-radiation  is  from  the  denser,  cooler  plasma. 

C.  Double-structure  Density  Gradient 

An  easy  way  to  simulate  numerically  a  prepulse  is  by  intro¬ 
ducing  a  double  density  gradient  scale  length.  We  chose  30  and  lOOp 

3 

for  the  standard  case.  Results  for  this  case  appear  in  Fig.  7  and  8. 

4 

Also  plotted  are  experimental  results  for  backscatter  for  laser  irra- 
diance  varying  between  10^  and  10^  W/cm^.  The  results  of  Ref.  4  are 
corrected  to  account  for  the  fact  that  our  simulations  are  with  a 
given  prepulse,  not  a  given  prepulse  to  main  pulse  ratio. ^ 

An  interesting  feature  of  these  results  is  the  maximum  elec¬ 
tron  temperature  observed  in  the  runs.  In  these  runs,  the  main  pulse 
is  of  a  short  pulse  type.  Once  backscatter  becomes  important,  the 
light  does  not  reach  the  critical  region  and  the  electron  temperature 
does  not  increase  any  more.  Also,  the  light  does  not  have  time  to 
burn  its  way  through  to  the  critical  density  region.  The  incident 
laser  energy  is  delivered  instead  to  the  ions  which  become  very 

g 

energetic  (2  x  10  cm/sec  or  100  keV  ions).  This  feature  could  be 
used  for  applications  which  will  be  mentioned  later. 

Another  dependence  of  the  double-structure  density  gradient 
was  investigated.  Instead  of  using  30p  on  the  density  gradient  scale 
length  near  the  critical  surface,  lOp  to  50p  was  tried  instead.  This 
variation  was  thought  to  mock-up  a  variation  in  strength  of  prepulse 
and/or  variation  in  delay  time.  Results  are  shown  in  Fig.  9.  The 
main  result  from  this  series  of  runs  is  the  small  variation  observed 
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due  to  changes  of  this  parameter.  This  feature  also  is  consistent 
with  the  experimental  results.  To  summarize,  both  our  simulations 
and  the  experiment  show  a  strong  increase  in  backscatter  with  increase 
in  irradiance  for  double  structure  pulses,  and  a  much  weaker  increase 
in  backscatter  as  strength  of  prepulse  to  main  pulse  increases  and/or 
delay  time  increases  (when  x  varies  from  30  to  50p  for  example). 

D.  Different  Elements 

Sensitivity  of  absorption  and  other  quantities  to  various 
elements  was  considered.  We  chose  aluminum  instead  of  polyethylene  as 
another  target  material  in  order  to  consider  the  largest  possible 
variation  in  target  element  without  introducing  unnecessary  compli¬ 
cations  due  to  the  radiation  processes.  It  was  expected  that  if  the 
temperature  reached  5  keV,  A1  should  be  fully  ionized  and  as  a  result, 
only  bremsstrahlung  losses  needed  to  be  taken  into  account,  like  in 

our  present  version  of  the  code.  A  test  was  made  with  A1  at  10^ 

2 

W/cm  for  a  short  pulse  and  a  double-structure  density  gradient.  The 
results  obtained  were  extremely  close  quantitatively  to  those  obtained 
for  a  CH2  target  for  all  the  quantities  plotted.  This  proves  that 
classical  inverse  bremsstrahlung  plays  a  very  insignificant  role  since 
it  is  five  times  stronger  at  that  laser  irradiance  for  aluminum  than 
for  CH2  but  the  absorption  and  temperature  stay  the  same. 


35 


* 


E.  Different  Laser  Wavelengths 

All  of  the  runs  reported  above  were  obtained  for  the  case  of 
Nd  laser  light  (X  =  1.06p).  Runs  were  also  performed  for  light  of 
other  wavelength.  Our  most  extensive  set  of  runs  is  for  a  CO, 

(X  =  10. 6u)  laser  produced  CH^  plasma.  To  do  so,  scale  lengths  and 
times  for  the  short  and  long  pulses  of  a  Nd  laser  plasma  were  multi¬ 
plied  by  ten  and  irradiances  were  divided  by  100.  Results  appear  in 
Fig.  10  and  11  for  short  pulses  and  in  Fig.  12  and  13  for  long  pulses. 
Overall,  results  are  quite  similar  to  those  obtained  for  Nd  light. 
Backscatter  is  observed,  both  in  the  cases  of  short  and  long  pulse 
experiments.  Comparable  levels  are  reached  for  slightly  larger  values 
of  IX2  for  CO2  than  for  Nd  laser  light,  and  for  the  short  pulse 
configuration,  absorption  by  ion  acoustic  turbulence  is  somewhat  less 

than  it  is  for  a  Nd  laser  plasma.  Notice  that  for  long  pulses  our 

13  2 

calculation  predicts  good  absorption  to  irradiances  beyond  10  W/cm  . 

To  test  the  effect  of  a  shorter  wavelength,  two  other  simu¬ 
lations  were  done,  one  at  0.53p  for  a  Cf^  plasma  at  an  irradiance  of 

15  2  52 

10  W/cm  which  we  could  compare  to  an  experiment,  and  one  at  0.26u 

16  2 

for  a  plasma  at  10  W/cm  .  In  the  first  case,  the  laser  pulse 
width  at  half-maximum  was  taken  to  be  70  psec  whereas  it  was  shortened 
to  40  psec  in  the  second  case.  The  initial  density  gradient  scale 
length  was  chosen  to  be  5u  but  it  was  observed  that  in  both  cases,  by 
the  end  of  the  pulse  it  was  more  like  15-20u.  For  the  0.53u  case, 
total  absorption  was  47%  vs.  45%  reported  experimentally.  The  maximum 
electron  temperature  obtained  numerically  was  2.9  keV  vs.  2  keV 
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measured  from  a  hard  X-ray  spectrum.  The  maximum  blowoff  velocity  was 

Q 

1.1  10  cm/sec  and  the  flux  limit  was  0.1.  It  was  also  found  in  our 
computation  that  a  fraction  of  23%  was  resonantly  absorbed  and  6%  only 
was  backscattered.  For  the  0.26m  case,  the  total  absorption  was 
further  enhanced  to  54%  with  the  contribution  of  the  resonant  absorp¬ 
tion  remaining  steady  at  23%.  The  specular  reflection  reached  42%  and 
backscatter  diminished  to  4%.  The  maximum  temperature  was  6.5  keV, 

g 

blowoff  velocity  1.5  10  cm/sec  and  flux  limit  f  =  0.1.  Thus  as  the 
wavelength  is  shortened,  most  of  the  absorption  increase  comes  from  an 
enhancement  of  the  absorption  due  to  ion  acoustic  instability  and 
inverse  bremsstrahlung. 

As  a  final  example  we  consider  millimeter  sources.  For 

instance  the  NRL  gyrotron  produces  about  150  kW  at  a  wavelength  of 
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8mm.  The  critical  density  is  n£  =  1.6  x  10  cm  .  Assuming  that 

the  beam  can  be  collimated  to  a  5  cm  radius,  this  corresponds  to  an 

irradiance  of  2  x  10  W/cm  ,  which  is  about  20  times  higher  than  the 

f  =  1  thermal  flux  n  mV  3  for  a  1  eV  plasma.  We  have  simulated  this 

c  e 

case  for  a  1  eV  hydrogen  plasma  for  a  'short  pulse',  where  the  gradient 
scale  length  is  8  cm  and  the  pulse  width  is  600  nanoseconds.  The 
total  absorption  is  35%,  resonant  absorption  is  29%,  specular  reflec¬ 
tion  is  65%  and  there  is  no  backscatter.  At  the  critical  density  the 
electron  temperature  increased  to  about  18  eV  but  there  was  no  heating 
anywhere  else.  The  return  current  driven  ion  acoustic  instability  was 
not  excited  anywhere.  Although  it  is  unphysical,  we  also  simulated  a 
'long  pulse'  at  2  x  10  W/cm  with  an  initial  density  gradient  scale 


length  of  80  cm  for  a  total  time  of  0.8psec.  (This  is  unphysical 
because  the  beam  radius  is  less  than  one  tenth  of  the  gradient  scale 
length;  two  dimensional  effects  should  be  very  important).  The 
absorption  is  much  greater,  about  79%,  but  it  is  all  resonant  absorp¬ 
tion  (17%)  and  inverse  bremsstrahlung;  again  the  return  current  driven 
ion  acoustic  instability  is  not  excited  anywhere.  Thus  our  simulations 
give  the  result  that  millimeter  wave  plasma  interactions  (and  presum¬ 
ably  microwave  plasma  interactions  also)  do  not  give  rise  to  a  return 
current  driven  ion  acoustic  instability. 

VI.  Suggestions  for  Further  Experiments 

A.  Varying  Pulse  Shape  to  Optimize  Absorption 

We  have  attempted  to  optimize  the  light  absorption  by  varying 

the  laser  pulse  shape.  In  all  cases  the  laser  system  was  operated  in 

a  long  pulse  mode  with  an  initial  density  gradient  scale  length  of 

15  2 

100pm  and  an  average  ir radiance  of  1.5  x  10  W/cm  .  Two  periods  of 
the  six  pulse  shapes  we  used  are  shown  in  Fig.  14.  For  a  uniform  laser 
pulse  (a),  the  fractional  absorption  was  67%  after  the  system  reached 
steady  state. 

Next  we  studied  the  effect  of  varying  the  laser  irradiance 
with  time  by  selecting  several  periodic  pulses  with  period  of  lOOpsec. 
All  absorption  efficiencies  we  quote  are  averaged  over  a  time  of 
lOOpsec  which  was  taken  to  begin  after  the  initial  transients  died 
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away.  First  we  tried  to  see  whether  relatively  small  changes  in 

pulse  shape  could  have  a  significant  effect.  To  test  this  out  we 

tried  two  sawtooth  oscillations,  first  (B)  an  abrupt  rise  to 

2.25  x  10^  W/cm2  and  a  linear  fall-off  to  0.75  x  10^"*  W/cm2  in  a 

hundred  psec;  and  second,  (C)  an  abrupt  fall  to  0.75  x  10^  W/cm2 

15  /  o 

followed  by  a  linear  rise  to  2.25  x  10  W/cm2  followed  by  a 
linear  rise  to  2.25  x  10^  W/cm2  in  100  psec.  The  absorption 
efficiency  for  the  former  was  62%  and  for  the  latter  was  59%. 

Since  relatively  small  perturbations  in  the  waveform  did 
not  have  a  significant  effect  on  absorption,  we  tried  larger  pertur- 

14  ~ 

bations.  First,  the  irradiance  (D)  was  set  at  3  x  10  W/cm2  for 

the  first  20  psec  and  then  raised  to  1.8  x  10^“*  W/cm2  for  the  last 

80  psec.  Second,  the  irradiance  (E)  was  raised  to  5  x  10^  W/cm2 

14  , 

for  20  psec  and  then  lowered  to  .625  x  10  W/cm2  for  the  remaining 

80  psec.  Finally,  the  irradiance  (F)  was  raised  to  5  x  10^  W/cm2 

for  5  psec  and  then  lowered  to  (~)  x  10  W/cm2  for  the  remaining 

95  psec.  The  fractional  absorptions  were  respectively  44%,  46%  and 

56%.  Thus  violently  changing  the  pulse  shape  appears  to  be  quite 

harmful  for  fractional  absorption.  The  best  absorption  is  then 

predicted  to  result  from  as  smooth  and  clean  a  pulse  as  possible. 

This  result  may  have  application  to  recent  experiments 

which  tested  the  fractional  absorption  as  a  function  of  position  of 

6  7 

the  target  with  respect  to  the  focal  point  of  the  lens.  ’  The  best 
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absorption  occurred  for  the  target  in  the  focal  plane.  In  the  NRL 
experiment, ^  the  beam  quality  was  also  measured  as  a  function  of 
distance  from  focal  plane.  At  best  focus,  the  beam  quality  was  by  far 
the  best.  The  curves  of  constant  intensity  were  simple  closed  curves, 
one  closed  curve  for  each  intensity.  Further  from  the  focal  plane, 
the  laser  light  had  many  more  hot  spots  and  the  beam  quality  was  much 
worse.  If  these  hot  spots  in  space  behave  the  same  as  the  hot  spots 
in  time  which  we  simulated,  the  fractional  absorption  should  maximize 
for  the  target  at  the  focal  plane. 

B.  Further  Studies  of  Resonant  Absorption 

54 

Recently,  Manes  et  al  have  performed  short  pulse  experi¬ 
ments  to  study  resonant  absorption.  By  tilting  the  target  for  s  and 
p  polarized  light,  they  found  a  maximum  absorption  where  the  theory 
of  resonant  absorption  predicts.  According  to  our  simulations, 
resonant  absorption  is  the  dominant  mechanism  for  short  pulse  experi¬ 
ments.  However  for  long  pulse  and  double  structure  pulse  experiments, 
it  plays  a  less  important  role.  Therefore  we  feel  that  it  would  be 
interesting  to  repeat  this  experiment  where  the  laser  is  run  in  a 
long  pulse  or  double  structured  pulse  mode  at  an  irradiance  of  perhaps 
10^~>  W/cm2.  If  our  model  is  correct,  there  should  be  less  evidence 
for  resonant  absorption. 

C.  Production  of  a  Dense,  Hot  Ion  Plasma 

As  we  have  seen  in  the  last  section,  once  the  irradiance 
increases  above  10^  W/cm2  for  a  double  structured  laser  pulse,  the 


absorption  is  predicted  to  drop  drastically  due  to  the  Brillouin 


backscatter.  However  a  fraction  of  the  incident  laser  energy  equal  to 

roughly  2  ((T  +  3  T  )/Mc2)1/2  is  deposited  directly  into  the  ions  by 

the  backscatter  instability.  For  irradiances  of  order  1017  W/cra2  our 

code  predicts  almost  complete  backscatter  and  an  ion  temperature  of 

ten  to  twenty  keV  for  the  entire  underdense  plasma.  (However  the 

electron  temperature  is  much  smaller). 

This  result  suggests  that  it  might  be  possible  to  produce  a 

dense  (ng  'v  10^  -  10^  cm  3)  plasma  with  very  high  ion  temperature. 

We  do  not  believe  any  other  plasma  source  at  this  density  and  ion 

temperature  exists.  For  one  thing  this  plasma  could  serve  as  a  very 

compact  thermal  source  of  thermonuclear  reaction  products.  To  test 

this  idea,  we  have  run  our  code  on  a  double  structured  plasma  pro- 

20 

duced  by  a  prepulse.  The  scale  length  is  100  pm  for  n  <10  and  is 

e 
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30pm  for  10  <  ng  <  10  .  This  plasma  is  then  irradiated  with  a 

75  psec  FWHM  Gaussian  pulse  with  a  maximum  irradiance  of  2  x  10^7 

W/cm2,  like  that  discussed  in  Section  V.  For  a  D-T  plasma,  the  ion 
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temperature  reaches  a  maximum  of  15  keV  and  for  n^  <  10  cm  is 

everywhere  greater  than  7keV  at  the  peak  of  the  pulse.  The  total 
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number  of  neutrons  produced  is  4.8  x  10  / cm‘  and  the  total  energy  in 

reaction  products  is  136  J/cm2  as  compared  to  laser  input  energy  of 
1.07  x  10^  J/cm2.  The  plasma  is  so  hot  and  dense  that  it  can  even 

serve  as  a  source  of  advanced  reaction  products.  For  instance  we  have 

,  12 
rerun  the  calculation  for  a  D-HeJ  plasma.  The  output  is  now  1.4  x  10 

protons/cm2,  or  3.94  J/cm2  of  reaction  products. 


Although  this  theoretical  result  is  admittedly  quite  specu¬ 
lative,  there  are  at  least  two  reasons  (other  than  our  code  results) 
to  believe  it  may  be  true.  First  of  all,  it  is  well  established 
theoretically  that  Brillouin  scattering  does  ultimately  deposit  energy 

2  V  /c  times  the  laser  energy  directly  into  the  ions.  Secondly,  the 
s 

4 

experimental  results  of  Ripin  and  McLean  show  that  for  a  double 

structured  plasma,  the  backscatter  increases  with  irradiance.  For 

instance  the  plot  of  backscattered  energy  versus  laser  energy  (Fig.  3 

of  Ref.  4)  shows  no  indication  of  saturation,  and  at  an  irradiance  of 
16  2 

10  W/cm  (the  highest  irradiance  shown)  50%  of  the  energy  is  back- 
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scattered.  If  this  trend  continues  for  irradiance  of  10  W/cm  , 
which  seems  likely,  and  if  the  backscatter  is  due  to  stimulated 
Brillouin  backscatter,  which  also  seems  likely,  then  there  is  really 
no  way  to  avoid  producing  this  dense  hot  ion  plasma.  On  the  other 
hand,  for  such  a  large  laser  energy  input,  the  value  of  the  trapping 
width  compared  to  thermal  velocity  (Eq.  (8))  was  calculated  to  be 
much  larger  than  for  our  other  simulations.  Hence  the  ion  wave  may 
be  driven  into  a  highly  nonlinear  regime,  and  conventional  wave¬ 
breaking  could  also  be  playing  an  important  role.  In  any  case,  it 

seems  to  us  that  irradiating  a  prepulse  produced  plasma  with  a  main 
17  2 

pulse  of  10  W/cm  would  be  a  very  interesting  experiment. 

VII.  Conclusions 

We  conclude,  from  the  evidence  presented  here,  that  fluid 
simulation  with  anomalous  transport,  is  an  accurate  and  cost  effec¬ 
tive  way  to  describe  the  absorption,  backscatter  and  flux  limitation 


in  a  laser  produced  plasma.  Our  calculations  give  quantitatively 
correct  absorption  of  laser  light  over  five  orders  of  magnitude  in 
irradiance,  for  short  pulses,  long  pulses,  and  double  structured 
pulses. 

There  are  at  least  two  results  of  our  calculation  which  may  have 
direct  application  to  laser  fusion.  First,  the  absorption  appears  to 
maximize  for  smooth,  clean  pulses.  This  obviously  indicates  that  a 
high  quality  laser  and  focusing  system  will  give  best  results.  Second, 

if  low  irradiance,  long  pulses  prove  to  be  a  viable  approach  to  laser 

5  55  13  2 

fusion,  *  a  CC^  laser  with  an  irradiance  of  about  10  W/cm  could 

be  very  attractive.  Although  the  absorption  efficiency  is  only  about 

half  that  of  a  Nd  laser  at  the  same  irradiance,  its  efficiency  is 

much  more  than  twice  as  great,  so  more  of  the  input  energy  to  the 

laser  is  coupled  to  the  pellet.  The  maximum  electron  temperature  for 

this  CO2  laser  produced  plasma  is  about  5  keV.  However  this  is  the 

21  -3 

temperature  only  at  the  critical  density.  At  n  'v  10  cm  ,  the 
temperature  is  typically  around  1  keV. 

Acknowledgments 

Many  discussions  with  Dr.  B.  H.  Ripin  are  acknowledged  concern¬ 
ing  the  interpretation  of  experimental  data.  We  would  also  like  to 
thank  him  and  Dr.  F.  C.  Young  for  making  available  to  us  experimental 
data  before  publication  (Fig.  4).  We  would  like  to  thank  Drs.  S. 

Bodner  and  R.  Lehmberg  for  discussions  concerning  our  model  for  SBS. 
This  work  was  supported  by  the  U.  S.  Department  of  Energy. 


References 


1.  W.  M.  Manheimer,  D.  G.  Colombant  and  B.  H.  Ripin,  Phys.  Rev.  Lett. 
38,  1135  (1977). 

2.  W.  M.  Manheimer  and  D.  G.  Colombant,  Phys.  Fluids  21^,  1818  (1978). 

3.  B.  H.  Ripin,  F.  C.  Young,  J.  A.  Stamper,  C.  M.  Armstrong,  R. 
Decoste,  E.  A.  McLean  and  S.  E.  Bodner,  Phys.  Rev.  Lett.  39,  611 
(1977). 

4.  B.  H.  Ripin  and  E.  A.  McLean,  NRL  Memo  Report  3964  (1979);  to  be 
published  in  Appl.  Phys.  Lett. 

5.  B.  H.  Ripin,  R.  R.  Whitlock,  F.  C.  Young,  S.  P.  Obenschain,  E.  A. 
McLean  and  R.  Decoste,  NRL  Memo  Report  3965  (1979);  B.  H.  Ripin, 

R.  Decoste,  S.  P.  Obenschain,  S.  E.  Bodner,  E.  A.  McLean,  F.  C. 
Young,  R.  R.  Whitlock,  C.  M.  Armstrong,  J.  Grun,  J.  A.  Stamper, 

S.  H.  Gold,  D.  J.  Nagel,  R.  H.  Lehmberg,  and  J.  M.  McMahon,  to  be 
published  in  Phys.  Fluids. 

6.  B.  H.  Ripin,  NRL  Memo  Report  3684  (1977). 

7.  C.  G.  M.  Van  Kessel,  J.  N.  Olsen,  P.  Sachsenmaier ,  R.  Sigel,  K. 
Eidmann  and  R.  P.  Godwin,  presented  at  the  10th  European  Conf.  on 
Laser  Interaction  with  Matter,  Ecole  Polytechnique,  Palaiseau, 
France  (1976). 

8.  J.  P.  Freidberg,  R.  W.  Mitchell,  R.  L.  Morse  and  L.  J.  Rudinski, 
Phys.  Rev.  Lett.  28,  795  (1972). 

9.  D.  W.  Forslund,  J.  M.  Kindel,  K.  Lee,  E.  L.  Lindman  and  R.  L. 

Morse,  Phys.  Rev.  A  11,  679  (1975). 

10.  K.  G.  Estabrook,  E.  J.  Valeo  and  W.  L.  Kruer,  Phys.  Fluids  1£,  1151 
(1975). 


* 


44 


11.  D.  W.  Forslund,  J.  M.  Kindel  and  E.  L.  Lindman,  Phys.  Fluids  18, 
1017  (1975). 

12.  R.  A.  Haas,  W.  C.  Mead,  V.  L.  Kruer,  D.  W.  Phillion,  H.  N, 
Kornblum,  J.  D.  Lindl,  D.  MacQuigg,  V.  C.  Rupert,  and  K.  G. 
Tirsell ,  Phys.  Fluids  20,  322  (1977). 

13.  H.  D.  Shay,  R.  A.  Haas,  V.  L.  Kruer,  M.  J.  Boyle,  D.  W.  Phillion, 
V.  C.  Rupert,  H.  N.  Kornblum,  F.  Rainer,  V.  V.  Slivinsky,  L.  N. 
Koppel,  L.  Richards,  and  K.  G.  Tirsell,  Phys.  Fluids  21^,  1634 
(1978). 

14.  M.  D.  Rosen,  D.  W.  Phillion,  V.  C.  Rupert,  W.  C.  Mead,  W.  L. 
Kruer,  J.  J.  Thomson,  H.  N.  Kornblum,  V.  W.  Slivinsky,  G.  J. 
Caporaso,  M.  J.  Boyle  and  K.  G.  Tirsell,  Preprint  UCRL  -  82146 
(1978). 

15.  G.  Zimmerman  and  W.  L.  Kruer,  Comm.  Plasma  Phys.  2_,  85  (1975); 

G.  B.  Zimmerman,  Lawrence  Livermore  Laboratory  Report  No.  UCRL- 
74811,  1973  (unpublished). 

16.  W.  M.  Manheimer,  Phys.  Fluids  ^0,  265  (1977). 

17.  W.  M.  Manheimer,  C.  E.  Max,  and  J.  J.  Thomson,  Phys.  Fluids  21 , 

2009  (1978). 

18.  R.  Faehl  and  W.  L.  Kruer,  Phys.  Fluids  20_,  55  (1977). 

19.  A.  Mase  and  T.  Tsukishima,  Phys.  Fluids  _18,  464  (1975). 

20.  C.  C.  Daughney,  L.  S.  Holmes  and  J.  V.  M.  Paul,  Phys.  Rev.  Lett. 

25,  497  (1970). 

21.  K.  Muraoka,  E.  L.  Murray,  J.  M.  W.  Paul  and  D.  R.  Summers,  J. 
Plasma  Phys.  H),  135  (1973). 


45 


22.  Y.  Kowai  and  M.  Guyot,  Phys.  Rev.  Lett.  3j),  1141  (1977). 

23.  S.  M.  Hamburger  and  J.  Jancarik,  Phys.  Fluids  13,  825  (1972). 

24.  R.  E.  Slusher,  C.  M.  Surko,  D.  R.  Moler  and  M.  Porkolab,  Phys. 

Rev.  Lett.  _36,  674  (1976). 

25.  A.  Y.  Wong,  B.  H.  Quon  and  B.  H.  Ripin,  Phys.  Rev.  Lett.  30, 

1299  (1973). 

26.  B.  H.  Quon,  A.  Y.  Wong  and  B.  H.  Ripin,  Phys.  Rev.  Lett.  32, 

406  (1974). 

27.  M.  Lampe,  W.  M.  Manheimer,  J.  B.  McBride,  J.  H.  Orens,  R.  Shanny, 
and  R.  Sudan,  Phys.  Rev.  Lett.  26,  1221  (1971). 

28.  M.  Lampe,  W.  M.  Manheimer,  J.  B.  McBride,  J.  H.  Orens,  K. 
Papadopoulos,  R.  Shanny  and  R.  N.  Sudan,  Phys.  Fluids  15^,  662 
(1972). 

29.  M.  Lampe,  W.  M.  Manheimer,  J.  B.  McBride,  and  J.  H.  Orens,  Phys. 
Fluids  15,  2356  (1972). 

30.  M.  Lampe,  I.  Haber,  J.  H.  Orens  and  J.  P.  Boris,  Phys.  Fluids  17 , 
428  (1974). 

31.  J.  A.  Wesson,  A.  Sykes  and  H.  R.  Lewis,  Plasma  Phys.  L5»  49  (1973). 

32.  D.  Biskamp  and  R.  Chodura,  Phys.  Fluids  _16,  893  (1973). 

33.  D.  Biskamp,  K.  U.  von  Hagenow  and  H.  Welter,  Phys.  Lett.  39A,  351 
(1972). 

34.  C.  T.  Dum,  Phys.  Fluids  21,  956  (1978). 

35.  M.  Z.  Caponi  and  R.  C.  Davidson,  Phys.  Rev.  Lett.  31,  86  (1973). 

36.  A.  M.  Sleeper,  J.  Weinstock  and  B.  Bezzerides,  Phys.  Fluids  16, 

1508  (1973). 


46 


♦ 


37.  W.  M.  Manheimer  and  R.  Flynn,  Phys.  Fluids  _17,  409  (1974). 

38.  D.  Choi  and  W.  Horton,  Phys.  Fluids  17,  2048  (1974). 

39.  M.  N.  Rosenbluth,  Phys.  Rev.  Lett.  29^ ,  565  (1972). 

40.  K.  Nishikawa  and  C.  S.  Liu  in  Advances  in  Plasma  Physics,  Vol.  6 
Wiley  &  Sons  (1976). 

41.  M.  N.  Rosenbluth,  R.  P.  White,  and  C.  S.  Liu,  Phys.  Rev.  Lett. 

31,  1190  (1973). 

42.  W.  M.  Manheimer,  Phys.  Fluids  _17,  1634  (1974). 

43.  D.  W.  Phillion,  W.  L.  Kruer,  and  V.  C.  Rupert,  Phys.  Rev.  Lett. 
37,  1529  (1977). 

44.  C.  L.  Tang,  J.  Appl.  Phys.  38,  390  (1967). 

45.  W.  M.  Manheimer  and  H.  H.  Klein,  Phys.  Fluids  _17,  1889  (1974). 

46.  P.  M.  Campbell,  R.  R.  Johnson,  F.  J.  Mayer,  L.  V.  Powers,  and 

D.  C.  Slater,  Phys.  Rev.  Lett.  J39,  274  (1977);  R.  L.  Morse  and 

R.  C.  Malone,  Bull.  Am.  Phys.  Soc.  21,  1028  (1976). 

47.  J.  M.  Dawson  and  C.  Oberman,  Phys.  Fluids  j > ,  394  (1963). 

48.  K.  Lee,  D.  Forslund,  J.  M.  Kindel  and  E.  L.  Lindman,  Phys.  Fluids 
20,  51  (1977). 

49.  B.  H.  Ripin,  Appl.  Phys.  Lett.  j30,  136  (1977). 

50.  F.  C.  Young  and  B.  H.  Ripin,  Bull.  Am.  Phys.  Soc.  22,  1112  (1977) 

submitted  for  publication. 

51.  J.  P.  Anthes,  M.  A.  Palmer,  M.  A.  Gusinow  and  M.  K.  Matzen,  Bull. 
Am.  Phys.  Soc.  23 ,  777  (1978). 


47 

* 


52.  F.  Amiranoff,  R.  Benattar,  R.  Fabro,  E.  Fabre,  C.  Garban, 


C.  Popovics,  A.  Poqueruse,  R.  Sigel,  C.  Stenz,  J.  Virmont  and 
M.  Weinfeld,  in  Plasma  Physics  and  Controlled  Thermonuclear 
Fusion  (International  Atomic  Energy  Agency,  Vienna,  1978)  (to  be 
published) . 

53.  M.  Read,  R.  Gilgenbach,  R.  Lucey,  K.  R.  Chu  and  V.  L.  Granastein, 
to  be  published. 

54.  K.  R.  Manes,  V.  C.  Rupert,  J.  M.  Auerbach,  P.  H.  Y.  Lee,  and 
J.  E.  Swain,  Phys.  Rev.  Lett.  _39,  281  (1977). 

55.  NRL  Laser-Plasma  Interaction  Group,  NRL  Memo  Report  3890  (1978), 
ed.  B.  H.  Ripin. 


48 


♦ 


The  code  which  is  listed  in  the  next  appendix  is  a  time-dependent 
one  dimensional,  one  fluid  hydrodynamic  code  with  two  temperatures.  It 
usesan  eulerian  description  of:  the  fluid  equations  and  uses  FCT  as  its 
basis  for  the  convection  algorithm.  It  uses  time-splitting  techniques 
and  all  the  convective  terms  are  treated  explicitly  whereas  the 
diffusive  terms  are  treated  implicitly.  The  hydrodynamic  equations 
are  very  common  and  have  received  wide  attention  and  coverage  in 
the  literature.  Other  equations  are  not  so  common  and  although  they 
do  not  present  formidable  difficulties,  care  must  be  exercised  when 
solving  them  in  order  for  them  to  reflect  the  physics  as  accurately 
as  possible. 

In  this  code,  the  field  equations  received  a  special  treatment 
because  of  the  following  fact  as  mentioned  in  section  III  ;  the 
coupling  between  Brillouin  backscatter  and  the  incident  field  is  very 
strong  and  must  be  reflected  in  the  algorithm.  The  amplitude  of  the 
Brillouin  backscatter  depend  on  the  amplitude  of  the  incident  field 
intensity  which  is  being  solved  for  in  the  field  equations.  Fields 
at  a  previous  time  step  cannot  be  used  because  they  lead  to  an  unstable 
situation.  Where  the  incident  field  is  large,  the  Brillouin  back¬ 
scatter  is  large  and  as  a  result,  the  incident  field  drops  very 
rapidly.  On  the  next  time  step,  the  Brillouin  backscatter  is  very 
small  and  the  field  penetrates  much  more  readily.  One  more  time  step 
and  the  Brillouin  backscatter  is  large  again,  making  the  laser  light 


recede.  In  order  to  avoid  this  unphysical  type  of  behavior,  a 


marching  technique  was  used  for  solving  Eq.  (20)  where  the  coefficient 
for  Brillouin  backscatter  was  computed  at  the  previous  location  using 
the  field  just  calculated  there.  This  marching  technique  can  lead  to 
potential  problems  if  it  is  not  reversible.  A  marching  technique 
typically  uses  values  at  the  previous  spatial  point  to  advance  the 
solution  to  the  next  point.  Because  we  retrace  our  steps  backwards 
in  order  to  get  the  reflected  intensity  as  well,  the  marching  will 
take  place  using  values  ahead  of  the  solution  if  we  look  at  it  from 
the  point  of  view  of  the  incident  field  calculation.  This  leads 
easily  to  irreversible  solution  which  then  could  be  interpreted  as 
"numerical"  absorption  for  the  field  equations. 

The  solution  to  this  problem  for  the  backscatter  term  as  well  as 
for  the  absorption  term  was  to  use  the  same  value  when  going  forward 
as  when  coming  back.  These  terms  were  calculated  at  the  middle  of  a 
cell  and  whatever  the  direction  the  calculation  was  proceeding,  the 
change  in  the  term  in  that  cell  was  the  same.  Tests  on  reversibility 
were  performed  and  found  accurate  to  within  one  percent. 
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Fig.  1  -  Schematic  diagram  of  importance  of  inverse  bremsstrahlung  and 
ion  acoustic  turbulence  as  a  function  of  irradiance  (for  Nd  laser  light) 
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Fig.  7  -  Plot  of  absorption,  resonant  absorption,  specular  reflection, 
and  backscatter  as  a  function  of  irradiance  for  a  double  structured  pulse 
Nd  laser  plasma.  The  line  which  extends  from  lO*’’  -  10*°  W/cni  is  the 
NRL  data 
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Fig.  10  -  Plot  of  absorption,  resonant  absorption,  specular  reflection 
and  backscatter  as  a  function  of  irradiance  for  a  short  pulse  CO2  laser 
produced  plasma 
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Fig.  11  -  Plot  of  maximum  electron  temperature,  ion  blow-off  velocity 
and  flux  limit  at  n  =  n-  as  a  function  of  irradiance  for  the  plasma  of 
Fig.  10 
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Fig.  12  -  Plot  of  absorption,  resonant  absorption,  specular  reflection 
and  backscatter  as  a  function  of  irradiance  for  a  long  pulse  CO2  laser 
produced  plasma 
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Fig.  13  -  Plot  of  maximum  electron  temperature,  maximum  ion  blow-off 
velocity  (i-;  continues  to  accelerate  for  the  duration  of  the  simulation) 


and  flux  limit  at  n  =  n  as  a  function  of  irradiance  for  the  plasma  of 
Fig.  12  C 


Fig.  14  -  Picture  of  the  temporal  structure  of  the  six  pulse  shapes 

simulated 
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